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Abstract 

A proof of existence of stationary dark soliton solutions of the cubic-quintic non- 
linear Schrodinger equation with a periodic potential is given. It is based on the 
interpretation of the dark soliton as a heteroclinic on the Poincare map. 
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1 Introduction 



In the present paper we consider stationary solutions of the cubic-quintic non- 
linear Schrodinger equation (CQNLS) 



iipt + ^ + V{x)iP - 9l \iP\ 2 iP - = (1) 
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with an L-periodic symmetric potential: V(x) = V(x + L) = V(—x) and real 
constants g\ and g 2 . More specifically, we are interested in solutions which 
allow the representation ip(t,x) = e~ tuJt <f)(x) where the function solves the 
stationary equation 



<t>xx + V(x)<f> - 91 <P 3 - 92 <P 5 = (2) 

with V(x) = uj + V(x), and which can be identified as dark solitons due to 
their nonzero boundary conditions 



(f>(x) — > 4>±{x) as x — > oo (3) 

with the functions (f>±(x) being real, sign definite, and L-periodic solutions 
of (2). As it is clear, 4>(x) acquires a zero value at some point of space, and 
therefore without loss of generality can be searched real, what is taken into 
account in the passage from (1) to (2). 

The model (1) having a general character, it describes weakly dispersive and 
weakly nonlinear wave processes, recently attracted considerable attention in 
connection with its application to the mean-field theory of Bose-Einstein con- 
densates [1]. In this context, ip(x) is a macroscopic wave function, |-0(a;)| 2 de- 
scribed linear atomic density, and V(x) is an optical lattice created by standing 
laser beams. In particular, existence of spatially localized pulses (also referred 
to as bright solitons) has been recently addressed in Refs. [2,3]. In Ref. [3] 
families of the solutions were presented and significant differences in behav- 
ior of stationary modes of the standard cubic nonlinear Schrodinger equation 
(gi 7^ 0, (?2 = 0) and of the quintic nonlinear Schrodinger (QNLS) equation 
(<7i = 0, gi ^ 0), have been found. Dark solitons of the NLS equation with a 
periodic potential have also been discussed in the small amplitude limit [4] and 
for a general case it has been studied numerically in Ref. [5] (see also [6] and 
references therein) . The approach developed in [5] was based on the numerical 
study of Poincare map generated by Eq. (2) considered at instants nL, the 
dark solitons appearing as heteroclinics on the map. The aim of the present 
paper is to extend earlier studies, providing for the first time a rigorous proof 
of the existence of a dark soliton solution of Eq. (1). From the mathematical 
point of view, the strategy of proof combines in a novel way several techniques 
from the classical theory of ODE's (upper and lower solutions [8] and trunca- 
ture arguments) and dynamics of planar homeomorphisms (topological degree 
[9,10] and free homeomorphisms [11]). 

We will restrict the consideration to the case g 2 > only, what rules out any 
possibility of blowing up solutions. Then without loss of generality we can set 
g 2 — 1 through a rescalling, what is done in what follows. 
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2 Existence of one-signed periodic solutions 



We start with a proof of existence of one-signed periodic solutions. To this end 
we impose the condition of boundness of V(x): V m i n < V(x) < V max , consider 
uj > —Vraini and introduce the notations A 2 = uo + V m i n and \\ = uo + V max . 
As it is clear < A 2 < V(x) < A 2 .. Next we consider two stationary equations 
(J = l,2) 



<f>3,*x + - 9i4>) - 4>) = (4) 

Treating these equations as dynamical systems, one easily finds the (only) two 
nontrivial equilibria ±pj where 



These are the hyperbolic points +pj and — pj which are connected by the 
heteroclinic orbits, which explicit forms read 



p,-a ? - tanh(kjx) 

<t>j = i (6) 

\J p] + a] - p] tanh 2 (kjx) 



where 



a 3 = \ 2 P 2 j + o#l> alld k 3 = PjJ Pj + if (7) 



Let us call pi )2 the positive equilibria. As it is clear 0i(O) = 02(0) = and 
4>i{x) < (f>2(x) for x > 0. 

At this point some considerations are required about the general second order 
equation 



4>xx = f{x,4>) ( 8 ) 

with f(x, 0) continuous with respect to the both arguments and L-periodic in 
x. The following definition is classical (see for instance [8] and its references). 

Definition 1 A function a : [xq, +oo) — > K. such that a xx (x) > f(x,a) 
(a xx (x) < f(x,a)) for all x > x$ is called a lower fupperj solution of eq. 
(8). 
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Now we can formulate 



Proposition 1 p 1 and p 2 are respectively lower and upper solutions of equa- 
tion (2). Hence here exists an unstable L-periodic solution between them. 

Proof: Let us observe that for j = 1 



Pi,** + V(x) Pl - g lP \ - p\> \\ Pl - g lP \ -/of = (9) 

and similarly for j = 2. Hence, p\ < p 2 are a couple of well-ordered lower and 
upper solutions respectively, therefore there exists a periodic solution between 
them [8]. Such solution is unstable because the associated Brouwer index to 
the Poincare map is —1 (see for instance [10]). □ 

Therefore we have a positive L-periodic solution of Eq. (2), we designate it as 
<f>+(x), such that p\ < <f>+(x) < p 2 . By the symmetry of equation we also have 
a negative solution <f>-(x) = —<f>+(x). 

To give an example of a periodic solution, we consider the simplified QNLS 
model (1) with gi = and with the potential 



V(x) =p 4 [2- k 2 sn 2 (p 2 x,k)] 2 (10) 

where sn(x, k) is the Jacobi elliptic function, k G [0, 1] is the elliptic modulus, 
and p > (examples of the exact periodic solutions for the cubic nonlinear 
Schrodinger equation, g 2 = 0, and with a specific potential, were obtained 
in [14,15]). The respective positive definite solution reads 



4>+(x) = pdn(p 2 x, k) (11) 

and it corresponds to the frequency uj = p A {k 2 — 3). To verify the stability 
of in the sense of the dynamical system (2) we consider small deviation 
i]){x) = <f>(x) — 4>+{x) at x — > oo, whose dynamics in the leading order is 
governed by the equation 



iP xx -U(x,k)^ = (12) 

with 

U(x, k) = p 4 (4 - k 2 - 6k 2 sn(p 2 x, k) 2 + 4k i sn{p 2 x, k) 4 ) > 

> P\\ - k 2 ) > (13) 

Thus, the obtained function (f>+(x) is a hyperbolic periodic solution of Eq. (2). 
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Considering now ip + (x,t) = <f>+(x) exp(i(3 — k 2 )i) as a solution of Eq. (1), 
performing the stability analysis as in Ref. [14], only slightly modified due to 
presence of quintic nonlinearity, and taking into account that (j>+(x) > one 
verifies that ip + (x,t) is linearly stable in the sense of the evolution problem 

(!)• 

More sophisticated models allowing exact sign definite periodic solutions can 
be constructed using a kind of "inverse engineering" (i.e. by obtaining poten- 
tials starting with given periodic solutions) as it is explained in [6]. 



3 Existence of a dark soliton 



In this section we prove the existence of a heteroclinic orbit connecting the 
periodic solutions 0_ and (f> + . A battery of preparatory lemmas are necessary. 

Lemma 1 If <p : [xq, +oo) — > R is a bounded solution of eq. (8), then the 
derivative <p x is also bounded in [x ,+oo). 

Proof: By the hypothesis \(j>(x)\ < M, where M is a constant, for all x > x . 
Then, by the mean value theorem there exits x n G (nL, (n + 1)L) such that 
(f)((n+l)L)—(f)(nL) = <p x (x n )L. From here <p x (x n ) < 2M/Lforall n. Applying 
the mean value theorem one more time one obtains 



\<t>x{x) - <t>x{x n )\ < L max \f(x, (f>)\, \/x e (nL, (n + 1)L) (14) 

\<p\<~M 

because \(f) xx (x)\ < max^|< M \f(x, 0)|. Then 



\<j> x (x)\ < L max : \f(x,<f>)\ + Wx > x . (15) 

\4>\<M L 

□ 

The following lemma is a key ingredient in our main result. 

Lemma 2 [12] Let P : R 2 — > R 2 be an orientation preserving homeomorphism 
with a unique fixed point pl such that ^{1 — P,p L } ^ 1 . Then for any p G R 2 
one of the following possibilities holds 

i) P n (p ) — > Pl as n — > +oo 

ii) ||P n (j9 )|| oo as n — > +oo 
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Here, 7(1 — P,po} is the local index associated to the Brouwer degree of po as 
a fixed point of the homeomorphism P. The proof relies in a basic property of 
free homeomorphisms exposed in [11], namely the u;-limit set of a given point 
has to be a connected set of the fixed point set. 

Lemma 3 If f(x,y) is strictly increasing in y, there exists a most one L- 
periodic solution of (8). 

Proof By contradiction, let us assume that 2/1,2/2 are two different L-periodic 
solutions of (8). First, let us suppose that 2/1,2/2 intersect among themselves, 
that is, there should be ti,t 2 such that z(t) = yi(t) —y2(t) verifies z(t ) = = 
z(ti) and z(t) > for t G (ii^)- However, by substracting the corresponding 
equations and using that / is strictly increasing, we get that z should be 
convex in (ti,t 2 ), which is a contradiction. Therefore, 2/1,2/2 do not intersect 
and we assume without loss of generality that yi(t) > 2/2(0 for all t. Again, 
z should be convex in the whole real line, but this is impossible because it is 
periodic. □ 

With the help of these previous lemmas we are able to prove an abstract 
convergence result. 

Theorem 1 Let <p : [x , +00) — > R a bounded solution of (8). Let us assume 
that 



9f(x,y) , s 

mm > 16 

xe[o,L] Qy 

Then there exists an L-periodic solution (p(x) such that 

lim (|0(x) - <p(x)\ + \<p x {x) - <p x {x)\) = (17) 



Proof: Let us define m = infa;> Xo <f>(x) and M = sup x>XQ <f>(x), as well as the 
truncated function 



f(x,y) 



f(x,y), Vye[m,M] 

f(x,M) + f y (x,M)(y-M), Vy>M (18) 
f(x,m) + f y (x,M)(y-m), Vy < m 



f is strictly increasing in y. Note also that is a solution of the truncated 
equation 
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<l>xx = hx,<f>) (19) 

Obviously, ^va. y ^± OQ f(x, y) = ±00 uniformly in x. Hence, there exist constants 
a and f3, such that a < (5 and f(x,a) < < f(x,/3) for all x. Such a and (3 
is a well-ordered coupled L-periodic lower and upper solutions, so there exists 
an L-periodic solution of (19) with index —1. This solution is unique by the 
previous lemma. Then, by Lemma 2, <f>(x) must converge to ip(x) since Lemma 
1 excludes it). As <p(x) G [to, M], it is a solution of (8). □ 

Theorem 2 Let us consider bounded functions a, j3 : [xo, +00) — > K. verifying 

1) a(x) < (5{x), Vrr > a;o 

^ a xx (x) > f(x,a) and (3 xx (x) > f(x,f3), Vx > a; 
T/ien i/iere exzsfo a solution <f>(x) of (8) such that 

a(x) < <f)(x) < P(x) (20) 
If moreover, there exists x such that 

3) min df ^ V) > 

xe[0,L] Qy 
ye[inf x >;P a(x),Bup x > X0 0(x)] 

then there exists an L-periodic solution (p(x) such that 

lim (|0(x) - ^(x)| + \<t> x {x) - <p x (x)\) = (21) 

x— »+oo 

Besides, (p(x) is the unique L-periodic solution in the interval [ini x > Xo a(x), sup x>XQ f3(x)]. 

Proof. The first assertion is a classical result due to Opial [13]. The second 
conclusion is a corollary of Theorem 1. □ 

In order to apply the above results to our model (1), we observe that i) now 
f(x, y) = y 5 + giy 3 — V(x)y, ii) due to parity of the potential one can consider 
x > and extend the obtained solution (f>(x) as an odd function to x < 0, 
iii) The functions <fii,2(x) given by (6) satisfy the conditions 1) and 2) of 
the Theorem 2 where a(x) = 4>\{x) and j3{x) = (f>2(x). Hence, in order to 
prove that there exists a solution <f>(x) of (1) converging to <f>±(x), found in 
Proposition 1, as x — > ±00, one has to verify the condition 3) of Theorem 2. 

As x can be taken arbitrarily large, this last condition is equivalent to 

min \5y 4 + 3g 1 y 2 - V(x)\ > (22) 

xe[o,L] 
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Starting with the case g± > we observe that (22) is now equivalent to hp\ + 
3(?iPi — ^1 > 0- The straightforward analysis of this last inequality, which takes 
into account the link (5), the definition of Ai ; 2, the requirement uo > —V m in 
necessary for > 0, gives the following estimate for the frequency 

^ - 9i\Jgi + 4^ + 4V min > V max - 5V min (23) 

Thus (23) is a sufficient condition for the existence of dark solitons at non- 
negative g\. 

Considering now g\ < the constrain (22) is reduced to 5pf — 3|(?i|p| — ^2 > 
and subsequently to the following inequality constrain to the frequency 

8u + 2g\ + 10V mm - 2V max - 5 9l yfi^j + g\ + 4 V~ 

+3g 1 ^4u; + g 2 1 +W max > (24) 

which must be satisfied simultaneously with u > —V min . 

In both cases considered before, the conclusion is that there exists an explicitly 
computable u such that the CQNLS has a different dark soliton for any 

UJ > CO . 



4 A concluding example 

As a concluding remark we consider an example illustrating a dark soltion, 
as well as other concepts introduced in the paper. To this end we recall the 
potential (10) and construct a dark soliton which tends to + given by (11) 
(we thus consider now gi — 0). This cannot be done analytically, and that is 
why we employ numerics. An example is shown in Fig. 1. 

For numerical obtaining the dark soliton, we have use shooting method. To 
this end we observed that (12) is a Hill's equation and thus taking into account 
(13), from Floquet's theorem one we have that ip(x) = P(x) exp(— ax) where 
P{x) = P(x + 2K(k)/p 2 ) is a periodic function and K(fc) is a complete elliptic 
integral of the first kind. For given parameters k and p one can easily compute 
the respective Floquet's exponent a. In particular for our choice of k — 0.5 
and p — 1, a ~ 2.014298. Thus, starting with the points Xi n i = 2K(k)/p 2 
where n is an integer, where <f)(xi n i) = <f>+(0) —C and (f) x (xini) — —Ca, and by 
varying the parameter C one can meet the condition 0(0) = 0. An example of 
implementation of this procedure is shown in Fig. 1. 



8 



Fig. 1. The dark soliton - heteroclinic (solid line), lower and upper solutions pi t 2 
(dashed-dotted line), and the hyperbolic periodic solutions 4>±{x) given by (11), for 
the parameters k = 0.5 and p = 1 
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